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QQ ■ Abstract 

Building on earlier work, we have given in [29] a proof of existence of Abrikosov vortex 
(-H ' lattices in the Ginzburg-Landau model of superconductivity and shown that the triangular 

O |. lattice gives the lowest energy per lattice cell. After [25] was published, we realized that it 

proves a stronger result than was stated there. This result is recorded in the present paper. 
The proofs remain the same as in [29], apart from some streamlining. 
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1 Introduction 

Q\ , 1.1 The Ginzburg-Landau equations. The Ginzburg-Landau model of superconductivity 

describes a superconductor contained in 51 C K™, n — 2 or 3, in terms of a complex order 
parameter \P : Q — > C, and a magnetic potential A : O — » The key physical quantities for 

04 ' the model are 

• the density of superconducting pairs of electrons, n s :— \^>\ 2 ; 

• the magnetic field, B :— curl A; 

• and the current density, J :— \tti(^>Va^)- 

In the case n = 2, curlj4 :— OA 2 - — SA 1 - is a scalar- valued function. The covariant derivative V4 

' 0x1 0x2 ™ 

is defined to be V — iA. The Ginzburg-Landau theory specifies that the difference between the 
supperconducting and normal free energies in a state (^,^4) is 

En{9,A):= f \V A *\ 2 + \ curl A\ 2 + ^(1 - |*| 2 ) 2 , (1.1) 

where n is a positive constant that depends on the material properties of the superconductor. 
(In the problem we consider here it is appropriate to deal with Helmholtz free energy at a fixed 
average magnetic field b :— J n curl A, where |f2| is the area or volume of fi.) It follows from 
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the Sobolev inequalities that for bounded open sets ft, £q is well-dehned and C°° as a functional 
on the Sobolev space H . The critical points of this functional must satisfy the well-known 
Ginzburg-Landau equations inside ft: 

A A t> = k 2 (|^| 2 - 1)*, (1.2a) 



curl* curl A = Im(^V A ^)- (1.2b) 

Here = — V| and curl* are the adjoints of Va and curl. Explicitly, V4F = — divF + 
iA ■ F, and curl* F = curlF for n = 3 and curl* / = (J£-, — Jj^) f° r n = 2. 

It is well-known that there exists a critical value k c (in the units used here, n c — l/y/2), that 
separates superconductors into two classes with different properties: Type I superconductors, 
which have k < n c and exhibit first-order phase transitions from the non-superconducting state 
to the superconducting state, and Type II superconductors, which have n > k c and exhibit 
second-order phase transitions and the formation of vortex lattices. Existence of these vortex 
lattice solutions is the subject of the present paper. 

1.2 Abrikosov lattices. In 1957, Abrikosov p] discovered solutions of (|1.2|) in n = 2 whose 
physical characteristics n s , B, and J are periodic with respect to a two-dimensional lattice, while 
independent of the third dimension, and which have a single flux per lattice cel{^|. (In what 
follows we call such solutions, with n s and B non-constant, C-lattice solutions, or, if a lattice C 
is fixed, lattice solutions. In physics literature they are variously called mixed states, Abrikosov 
mixed states, Abrikosov vortex states.) Due to an error of calculation Abrikosov concluded 
that the lattice which gives the minimum average energy per lattice celjfl is the square lattice. 
Abrikosov's error was corrected by Kleiner, Roth, and Autler [21], who showed that it is in fact 
the triangular lattice which minimizes the energy. 

1.3 Results. In this paper we combine and extend the previous technique to give a complete 
and self-contained proof of the existence of Abrikosov lattice solutions. As in previous works, 
we consider only bulk superconductors filling all R 3 , with no variation along one direction, so 
that the problem is reduced to one on R 2 . To formulate our results, for a lattice C C R 2 , we 
denote by ft c and \ fl c \ the basic lattice cell and its area, respectively (for details see Section [3]). 
The flux quantization (see below) implies that 

1^1 ~, (1.3) 

where b is the average magnetic flux per lattice cell, b :— j^rrj f n r. curl A We note that due to 
the reflection symmetry of the problem we can assume that b > 0. We define 




«c(r):=W- I--57-T . (1-4) 



where /3(r) is the Abrikosov parameter, introduced in (|9.1j) below. We will prove the following 
results. 

Theorem 1.1. Let \b— n 2 \ <C 1 and (n — k c (t))(k 2 — b) > 0. Then for every lattice £ satisfying 
(jl.3[) with n = 1, the following holds 



2 Such solutions correspond to cylindrical samples. In 2003, Abrikosov received the Nobel Prize for this 
discovery 

3 Since for lattice solutions the energy over R 2 (the total energy) is infinite, one considers the average energy 
per lattice cell, i.e. energy per lattice cell divided by the area of the cell. 
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(I) The equations (|1.2[) have an C-lattice solution in a neighbourhood of the branch of normal 



(II) The above solution is unique, up to symmetry, in a neighbourhood of the normal branch. 
(Ill) For (k — k c (t))(k 2 — b) ^ 0, the solution above is real analytic in b in a neighbourhood of 



(IV) For k 2 > 1/2, the lattice shape for which the average energy per lattice cell is minimized 
approaches the triangular lattice as b — > k 2 , in the sense that the shape parameter, tc, of 
C (see Subsection 3.3 below) approaches T tr ianguiar = e I7r / 3 in C. 

Since their discovery, Abrikosov lattice solutions have been studied in numerous experimental 
and theoretical works. Of more mathematical studies, we mention the articles of Eilenberger 
|15) . Lasher [22], Chapman [10] and Ovchinnikov [25] . 

The rigorous investigation of Abrikosov solutions began soon after their discovery. Odeh [24] 
sketched a proof of existence for various lattices using variational and bifurcation techniques. 
Barany, Golubitsky, and Turski [8 applied equivariant bifurcation theory and filled in a number 
of details, and Takac [55] has adapted these results to study the zeros of the bifurcating solutions. 
Further details and new results, in both, variational and bifurcation, approaches, were provided 
by [13]. In particular, [14] proved partial results on the relation between the bifurcation 
parameter and the average magnetic field b (left open by previous works) and on the relation 
between the Ginzburg-Landau energy and the Abrikosov function, and [13] (see also [M]) found 
boundaries between superconducting, normal and mixed phases. 

Among related results, a relation of the Ginzburg-Landau minimization problem, for a fixed, 
finite domain and external magnetic field, in the regime of k — > oo, to the Abrikosov lattice 
variational problem was obtained in [3] [3] . 

All the rigorous results above deal with Abrikosov lattices with one quantum of magnetic 
flux per lattice cell. Partial results for higher magnetic fluxes were proven in [TU] [5] . 

After introducing general properties of (|1.2[) in Sections [2131 we prove the above theorem in 
Sections [MS 

Acknowledgements The second author is grateful to Yuri Ovchinnikov for many fruitful dis- 
cussions. A part of this work was done during I.M.S.'s stay at the IAS, Princeton. 

2 Properties of the Ginzburg-Landau equations 

2.1 Symmetries. The Ginzburg-Landau equations exhibit a number of symmetries, that is, 
transformations which map solutions to solutions: 
The gauge symmetry, 



solutions. 



K 



2 



(*(x),A(x)) i y {e vq(x) ^{x),A{x) + Vr?(x)), V?y G C 2 (M 2 ,M); 



(2.1) 



The translation symmetry, 



(*(x), A{x)) i ^ (*(.t + t),A{x + t)), 



Vi G R 2 ; 



(2.2) 



The rotation and reflection symmetry, 



{V{x),A(x)) ^ (*(R- 1 x),RA(R- 1 x)) 



Vi?. g 0(2). 



(2.3) 
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2.2 Flux quantization. One can show that under certain boundary conditions (e.g., 'gauge- 
periodic', see below, or if SI = R 2 and £q < oo) the magnetic flux through Q is quantized: 

:= / curlA = 27m (2.4) 
Jn 

for some integer n. 

2.3 Elementary solutions. There are two immediate solutions to the Ginzburg-Landau 
equations that are homogeneous in . These are the perfect superconductor solution where 
\& g = 1 and Ag = 0, and the normal (or non-superconducting) solution where "f^r = and 
An is such that cur\A^ =: b is constant. (We see that the perfect superconductor is a solution 
only when &(A) — 0. On the other hand, there is a normal solution, (^/y = 0, An, curlAjv = 
constant), for any condition on $(/!).) 

Moreover, for any integer n there is a (n— ) vortex solution of the form 

*W(ar) = f n (r)e m0 and (a) = a„ (r)V(n0) , (2.5) 

where (r, 0) are the polar coordinates of i € I 2 , unique up to symmetry transformations (see 
[SOU). Note that = n . 



3 Lattice states 

3.1 Periodicity. Our focus in this paper is on states ("J, A) defined on all of R 2 , but whose 
physical properties, the density of superconducting pairs of electrons, n s :— | \I/ 1 2 , the magnetic 
field, B := curl^4, and the current density, J :— Im( , ^ , Y4 , i , ) , are doubly-periodic with respect to 
some lattice C We call such states C— lattice states. 

One can show that a state (*, A) e Hl oc (R 2 ; C) x H^ C (R 2 ; R 2 ) is a ^-lattice state if and only 
if translation by an element of the lattice results in a gauge transformation of the state, that is, 
for each t € £, there exists a function gt £ Hf or (R 2 ; R) such that 

*(a; + 1) = e iat{x) ^{x) and A(x + t)= A(x) + Vg t (x) 

almost everywhere. 

It is clear that the gauge, translation, and rotation symmetries of the Ginzburg-Landau 
equations map lattice states to lattice states. In the case of the gauge and translation symmetries, 
the lattice with respect to which the solution is periodic does not change, whereas with the 
rotation symmetry, the lattice is rotated as well. It is a simple calculation to verify that the 
magnetic flux per cell of solutions is also preserved under the action of these symmetries. 

Note that (^,^4) is defined by its restriction to a single cell and can be reconstructed from 
this restriction by lattice translations. 

3.2 Flux quantization. The important property of lattice states is that the magnetic flux 
through a lattice cell is quantized, i.e. (12.4[) holds, with il any fundamental cell of the lattice. 

Indeed, if \^>\ > on the boundary of the cell, we can write \& = |\I/|e and < 8 < 2tt. 
The periodicity of n s and J ensure the periodicity of V6* — A and therefore by Green's theorem, 
J Q curl A = <f gQ A = <f an V8 and this function is equal to 27m since \P is single- valued. 

Equation (|2.4[) then imposes a condition on the area of a cell, namely, (II. 3[) . 

3.3 Lattice shape. In order to define the shape of a lattice, we identify x € R 2 with 
z — x\ + ix2 <E C, and view £ as a subset of C. It is a well-known fact (see [4j) that any lattice 
£CC can be given a basis r, r' such that the ratio t ~ — satisfies the inequalities: 



(i) M > 1; 
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(ii) Imr > 0; 

(iii) -\ < Rer < \, and Rer > if |r| = 1. 

Although the basis is not unique, the value of r is, and we will use that as a measure of the 
shape of the lattice. 

Using the rotation symmetry we can assume that C has as a basis {rei,rr}, where r is a 
positive real number and e\ — (1, 0). 

4 Fixing the gauge and rescaling 

In this section we fix the gauge for solutions, (^/, A), of (|1.2|) and then rescale them to eliminate 
the dependence of the size of the lattice on b. Our space will then depend only on the number 
of quanta of flux and the shape of the lattice. 

4.1 Fixing the gauge. The gauge symmetry allows one to fix solutions to be of a desired 
form. Let Aq(x) = %Jx, where J is the symplectic matrix 



We will use the following preposition, first used by [24] and proved in [28] (we provide an alternate 
proof in Appendix |A| . 

Proposition 4.1. Let ("J, A) be an C-lattice state, and let b be the average magnetic flux per 
cell. Then there is a C-lattice state (<f>, Aq+o), that is gauge- equivalent to a translation of(^,A), 
such that 

(i) (f>(x + t) = e^ x ' Jt cf>(x) and a(x + 1) = a(x) for all t in a fixed basis of C; 

(ii) a has mean zero: J^a = 0; 

(iii) a is divergence-free: diva = 0. 

4.2 Rescaling. Suppose, that we have a /^-lattice state (^>,A), where C has shape r. Now 
let b be the average magnetic flux per cell of the state and n the quanta of flux per cell. From 
the quantization of the flux, we know that 



Let L T be the lattice spanned by r T and r T r, with Q T being a primitive cell of that lattice. We 
note that |f2 T | = 2irn. We summarize the effects of the rescaling above: 





(4.1) 



We set a :— 



(f ) 2 = -pF- We now define the rescaled fields (ip, a) to be 



(ip(x),a(x)) := (a 1 i'(ax),o-A(ax)). 



(A) (ip, a) is a £ T -lattice state. 

(B) ^Enr^A) = £x(i>,a), where A = and 

= K^TT2 / (Wl 2 + |curlo| 




(4.2) 
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(C) \P and A solve the Ginzburg-Landau equations if and only if ip and a solve 



(-A a - A)</> 



(4.3a) 



curl* curl a = Im^Vji/j) 



(4.3b) 



,2 



for A = The latter equations are valid on f2 r with the boundary conditions given in 
the next statement. 

(D) If (^,A) is of the form described in Proposition ^. 1[ then ip and a satisfy 

(a) i/j(x + t) = e^ x ' Jt ip(x) and a{x + t) = a{x), for t = r r , r T r, where a — Aft + 



(c) diva = 0. 

Our problem then is, for each n = 1, 2, . . ., find (ip, a), solving the rescaled Ginzburg-Landau 
equations (|4.3[) and satisfying (iv), and among these find the one that minimizes the average 
energy £ x . 

5 Asymptotics of solutions to (14.31) 

In this section we derive properties of families of solutions of (|4.3j) depending on b, in the regime 
of b — > n 2 n, provided such families exist. It might be convenient to reparametrize such families 
assuming that b, or A, depends on a parameter e — > 0. Most of the results of this section were 
first stated in [1] (see also [10]). In what follows we use the notation 



for the average of a function / over the lattice cell f2 T .The main result of this section is the 
following 

Proposition 5.1. Assume the equations (|4.3j) have a family, (ip e , a e , A c ), e — > 0, of Abrikosov 
lattice solutions, with (curla c ) = n, satisfying 



(with the first and second derivatives of the remainders obeying similar estimates). Then ij)®, a\ 
and Ai satisfy the equations 



a, with A$(x) := f 
(b) far a = 0, 



Jx, 




^ = «/>o + 0(e 3 ), a e = Aq + e 2 ai + 0(e 4 ), A e = n + e 2 Ai + 0(e 4 ), 



(5.1) 



A-A™ipo — and curlai 




(5.2) 



and 




(5.3) 



where /3(ipo) is the Abrikosov parameter, defined by 



KM ■ 



(IV'ol 4 ) 
(IV'ol 2 ) 2 - 



(5.4) 
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Furthermore, the energy is expressed as 



£a« Waffle) = IT + 



(« 2 - |)/3(^ ) + | 



{\Mr + 0(e b ). (5.5) 



Proof. Plugging (15.11) into (|4.3[) and taking e — > gives the first equation in (15.21) and 

curl* curl ai = Im('0oV J 4«V ; o)- (5-6) 
We show now that 

Im(^ V A ^o) - ~\ curl* |^o| 2 - (5.7) 

(Recall, that for a scalar function, f(x) £ R, curl* / = (c^/, —dif) is a vector.) It follows from 
(|6.3[) , Section [51 that Vo satisfies the first order equation 

((Va S )i+*(Va S ) 2 )Vo = 0. (5.8) 

Multiplying this relation by ^q, we obtain ^o(^A")i^o + *V'o(Va")2'0o = 0. Taking imaginary 
and real parts of this equation gives 

Im^o(VAj)iV'o = -Re-0 o (V j4 ^)2V'o = -^JVol 2 

and 

1 



2 



which, in turn, gives (|5.7J) . 

The equations (|5.6I) and (|5.7p give curlai = H — \\i^o\ 2 : with iJ a constant of integration. 
H has to be chosen so that J nT curia! = 0, which gives the second equation in (15. 2[) . 

Now we prove (|5.3p . We multiply the equation (|4.3a[) scalar ly (in L 2 (Vt T )) by tpo, use that 
the operator — A a is self-adjoint and (— A a — njipo — 0, substitute the expansions (|5.1I) and take 
e = 0. to obtain 



Ai / hfor + 2i/ ^ ai • Va»Vo + ^ / |^o|=0. (5.9) 



This expression implies that the imaginary part of the second term on the left hand side of (15. 9p 
is zero. (We arrive at the same conclusion by integrating by parts and using that divai = 0.) 
Therefore 

2i / ^Qdi • \7^n^ = -2 / ai • Im^oV^Vo) = -2 / ai • curl* curlai. 
Jo t ./n T ./n T 

Integrating the last term by parts, we obtain 2i J fiT ^o a i ' V4jV'o = ~2 /^(curlai) 2 . Using this 
equation and the second equation in (|5.2j) . we obtain 



2i/ 4> a ai-V A nJo = -\ [ \M A +\{\i>o?) I I V>o 1 2 - (5.10) 
This equation together with (|5.9[) and the definition (|5.4p gives (|5.3p . 
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Now, we prove statement (|5.5|) about the Ginzburg-Landau energy. Multiplying (I4.3a|) 
scalarly by ip and integrating by parts gives 

\V a ^\ 2 = n 2 [ (AH 2 -« 2 M 4 ). 
Jar 

Substituting this into the expression for the energy, we find 

..2 



£x(sl>,a) = 'io(^~ vd^l 4 ) + (I curla| 2 )). (5.11) 



K , A K 

Using the expansions (|5.ip and the facts that curl Aft = n and (curlai) = gives 



£ x Me,a t ) = y + ^ + ^ (-y (|^o| 4 ) + (I curl ai | 2 )) + 0{e% (5.12) 



Next, using the second equation in ()5.2[) together with the last equation in (|5.1|) and the 
definition (pT3]l yields ([5\5]l . □ 



Eqn (|5.3p fixes the parameter e uniquely up to the normalization of ipo. Indeed, we observe 

A-71 

definition A = ^ and (J5T3J), yields 



that the third equation in (|5.1[) implies e 2 = ^r- 2 - + 0((A — n) 2 ), which, together with the 



* = w^Ml imP> +oll *- bn (5 ' 13) 

This equation implies the following necessary condition on existence of the solutions: 

b < k 2 if (n 2 - i)/3(Vo) + \ > and 6 > k 2 if (k 2 - i)^ ) + \ < 0. (5.14) 
(j57T3)) together with (j5~5]) yields 

ft. (*,„,) = ^ + 6- + 1 + <*(«• - »)•). (5.15, 

6 The linear problem 

In this section we solve the linear problem: — A^n^o = 'T'V'Oj for V' satisfying the gauge - periodic 
boundary condition ip(x + 1) = e~ x ' Jt ip(x), for t = r T , r T r, (see (|5.2[) and |D|).) We introduce 
the harmonic oscillator annihilation and creation operators, a n and (a™)*, with 

Tl ZTl 

a" := (Vi;)i + i{^A^2 = d Xl + id X2 + -x x + yx 2 . (6.1) 

One can verify that these operators satisfy the following relations: 

1. [a", (a™)*] = 2 curl Aft = 2n; 

2. -A^n - n= (a")*a n . 

As for the harmonic oscillator (see for example |18j). this gives explicit information about 

a{-A AS ) = {{2k+ l)n : fc = 0, 1, 2, . . .}, (6.2) 
and each eigenvalue is of the same multiplicity. Furthermore, the above properties imply 

nuU(-AA» -n) = nulla™. (6.3) 
We can now prove the following. 
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Proposition 6.1. null(— A^n — n) is given by 




mx\\{-l\ K -n) = {e^ X2(xi+lX2) ^ c k e ikV2 *^ {xi+lX2) \ c k+n = e ln * T e l2k7TT c k } (6.4) 

k— — oo 

and therefore, in particular, dime null(— Aaj — n) = n. 

Proof. We find null a n . A simple calculation gives the following operator equation 

e fi^i 2 a ™ e -f ki 2 =d Xl +id X2 . 

This immediately proves that tp £ nulla™ if and only if £ = e^^ip satisfies d Xl £ + id X2 £ = 0. 
We now identify x £ M 2 with z — x\ + ix% £ C and see that this means that £ is analytic. We 
therefore define the entire function to be 

0(z) = £ 

The quasiperiodicity of ■0 transfers to as follows 

6(Z + 7T) = 0(2!), 

0(z + ttt) = e- 2inz e-" l7rr 0(z). 

To complete the proof, we now need to show that the space of the analytic functions which 
satisfy these relations form a vector space of dimension n. It is easy to verify that the first 
relation ensures that have a absolutely convergent Fourier expansion of the form 

oo 

©(z)= c ^ 2klz - 

The second relation, on the other hand, leads to relation for the coefficients of the expansion. 
Namely, we have c k + n = e m7TT e 2kl ' KT c k and that means such functions are determined solely by 
the values of Co, . . . , c„_i and therefore form an n-dimcnsional vector space. □ 

7 Reformulation of the problem 

In this section we reduce two equations (|4.3p for tp and a to a single equation for ip. We 
introduce the spaces _S? n (r) := L 2 (fl T ,C) and if(r) := {a £ L 2 (fT,R 2 ) | (a) = 0, diva = 0, in 
the distributional sense}, with the domains being the Sobolev space of order two, J^(t), whose 
elements satisfy the quasiperiodic boundary condition (d) in part (|D]) above, and the Sobolev 
space of order two, J4?(t), whose elements satisfy the periodic boundary conditions with respect 
to C T , have mean zero, and are divergence free. Substituting a — Aq + a, we rewrite (|4.3[) as 

(£" - A)V> + 2ia ■ V A ^iP + |ct| 2 -0 + k 2 |^|V = 0, (7.1a) 

(M+\ 1 p\ 2 )a-lm(^V A ^) = 0, (7.1b) 

where 

L n := -A A « and M := curl* curl, (7.2) 

defined on the spaces S£ n (j) and -S?(t). Their properties that will be used below are summarized 
in the following propositions: 
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Proposition 7.1. L n is a self-adjoint operator on J4^(t) with spectrum cr{L n ) — { (2k + l)n : 
k = 0, 1, 2, . . .} and dime null(L" — n) = n. 

Proposition 7.2. M is a strictly positive operator on J$?(t) with discrete spectrum. 

The proofs of these results are standard and, for the convenience of the reader, are given 
below. 

Proof of Proposition \ 7. S\ The fact that M is positive follows immediately from its definition. 
We note that its being strictly positive is the result of restricting its domain to elements having 
mean zero. □ 



Proof of Proposition \ 7. 1\ First, we note that L n is clearly a positive self-adjoint operator. To see 
that it has discrete spectrum, we first note that the inclusion H 2 <-} L 2 is compact for bounded 
domains in M 2 with Lipschitz boundary (which certainly includes lattice cells). Then for any z 
in the resolvent set of L n , (L n — z)~ x : L 2 — > H 2 is bounded and therefore (L n — z)^ 1 : L 2 — » L 2 
is compact. In fact, the spectrum of L n was found explicitly in the previous section. This 
completes the proof of Proposition [73] □ 



We first solve the second equation (|7.1bl) for a in terms of ip, using the fact that M is a strictly 
positive operator, and that div J a = 0, (J a ) = 0, where J a :— Im{^^»+ a T/'}. The last two 
relations follow for any solution (ip, a) of (|4. 3|) by differentiating the equation £\(e sx ip, a+sWx) — 
£\(tp, a), w.r.to s at s — 0, which gives d^,E\(ip, a)ix^ + d a E\(ip, a)Vx = 0. Here d^£\(ip, a) and 
d a £\(ip,a) are the Gateaux derivatives of £\(ip,a) w.r.to ip and a. Since d^£\(ip,a) = 0, this 
yields = J QT (M — J a ) ■ V% = f QT div J a \. Since the last equation holds for any x € Hi(£l T , K), 
we conclude that div J a — 0. Choosing x = h- x, Vh € M 2 , in the equation = f QT (M— J a )-Vx, 
we find (J a ) — 0. Now, (|7. lb)) can be rewritten as a fixed point problem a = M~ 1 J a , which has 
a unique solution in Jif(r). The latter can be rewritten as a — a(ip), where 

a(V) = (M + H 2 )- 1 Im^^). (7.3) 

We collect the elementary properties of the map a in the following proposition, where we identify 
M' n {f) with a real Banach space using ip O ip := (Reip,Imip). 

Proposition 7.3. The unique solution, a(ip), of (|7.1b[) maps Jff n (r) to J4?(t) and has the 
following properties: 

(a) a(-) is analytic as a map between real Banach spaces. 

(b) a(0) = 0. 

(c) For any 8 € K, a(e lS ip) = a(ip). 

Proof. The only statement that does not follow immediately from the definition of a is (a). It 
is clear that Im^X^™^) is real-analytic as it is a polynomial in tp and V^, and their complex 
conjugates. We also note that (M — z)" 1 is complex-analytic in z on the resolvent set of M, 
and therefore, (M + IV-"! 2 ) -1 is analytic, (a) now follows. □ 

Now we substitute the expression (|7.3[) for a into (|7.1a[) to get a single equation 

F(A,^) = 0, (7.4) 

where the map F:Kx ^ n { T ) ~^ -^n( T ) i g defined as 

F(A,V) = (L n ' - X)<i/j + 2ia(ip) ■ Va»V> + |a(V>)lV + « 2 |^|V (7-5) 

For a map F(tp), we denote by d^F((f>) its Gateaux derivative in tp at <fr. The following 
proposition lists some properties of F. 
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Proposition 7.4. 

(a) F is analytic as a map between real Banach spaces, 

(b) for all X, F(X, 0) = 0, 

(c) for all X, d 4 ,F(X, 0) = L n - A, 

(d) for all S G R, F[X, e i& ip) = e lS F(X, tp). 

(e) for all ip, (ip,F(X,tp)) € R. 

Proof. The first property follows from the definition of F and the corresponding analyticity of 
a(ip). (b) through (d) are straightforward calculations. For (e), we calculate that 

(V, F(X, V0> = ty, (L n - A)V> + 2i f $a(j>) ■ VV> 



+ 2 f (a^)-A%M 2 +[ |a(^)| 2 |V| 2 + « 2 / NM 4 - 

The final three terms are clearly real and so is the first because L n — X is self-adjoint. For the 
second term we calculate the complex conjugate and see that 



2i / -)pa(ip) ■ Vip = -2i / ipa(tp) ■ V<ip = 2i (Vip ■ a(ij>))ij), 

JQJ JW JQJ 

where we have integrated by parts and used the fact that the boundary terms vanish due to 
the periodicity of the integrand and that div a(ip) — 0. Thus this term is also real and (e) is 
established. □ 



8 Reduction to a finite-dimensional problem 

In this section we reduce the problem of solving the equation .F(A, ip) = to a finite dimensional 
problem. We address the latter in the next section. We use the standard method of Lyapunov- 
Schmidt reduction. Let X := J? n (f) and Y := Jz?„(t) and let K = null(L" - n). We let P be 
the Riesz projection onto K , that is, 

P-=~ l{L n -z)- l dz, (8.1) 
2m / 7 

where 7 C C is a contour around n that contains no other points of the spectrum of L n . This 
is possible since n is an isolated eigenvalue of L n . P is a bounded, orthogonal projection, and if 
we let Z := nullP, then Y — K © Z. We also let Q :— I — P, and so Q is a projection onto Z. 
The equation F(X, if>) = is therefore equivalent to the pair of equations 

PF(X, Pi/j + Qi/>) = 0, (8.2) 
QF(A,PV + W) = 0. (8.3) 



We will now solve (I8.3[) for w = Qtp in terms of A and v = Pip. To do this, we introduce 
the map G:RxifxZ^Ztobe G(A, v, w) :— QF(X, v + w). Applying the Implicit Function 
Theorem to G, we obtain a function w : M x K — > Z , defined on a neighbourhood of (n, 0), such 
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that w = w(X,v) is a unique solution to G(X,v,w) = 0, for (A, v) in that neighbourhood. This 
solution has the following properties 

w(X,v) real-analytic in (X,v); (8.4) 



w(X,v) = 0(\v\ 2 ) and d x w(X,v) = 0{\v\ 2 ). (8.5) 

The last property follows from the fact that the last three terms in (|7.5[) are at least quadratic 
in ip = v + w. 

We substitute the solution w = w(X,v) into (|8.2p and see that the latter equation in a 
neighbourhood of (n, 0) is equivalent to the equation (the bifurcation equation) 

>y(X,v) := PF(X,v + w(X,v)) = 0. (8.6) 

Note that 7 : R x K — > C. We have shown that in a neighbourhood of (n, 0) in 1 x AT, (A, "0) 
solves F(X,tfj) — if and only if (A, v), with v — Pip, solves (|8.6[) . Moreover, the solution -0 of 
F(X, ip) = can be reconstructed from the solution v of (18. 6[) according to the formula 

tp = v + w(X,v). (8.7) 

Finally we note that w and 7 inherit the symmetry of the original equation: 

Lemma 8.1. For every i5 6 E, w(X, e lS v) — e lS w(X 7 v) and j(X, e l5 v) — e l<5 7(A, v). 

Proof. We first check that w(X, e l& v) = e l6 w(X, v). We note that by definition of w, 

G{X,e iS v,w(X,e lS v)) =0, 

but by the symmetry of F, we also have G(X, e zS v, e lS w{X, v)) — e lS G(X,v,w(X,v)) — 0. The 
uniqueness of w then implies that w(X, e lS v) = e lS w(X,v). We can now verify that 

7 (A, e lS v) = PF(X, e lS v + w(X, e iS v)) 

= e lS PF{X,v + w(X,v))) = e l<5 7 (A,f). □ 

Solving the bifurcation equation (|8.6p is a subtle problem, unless n — 1. In the latter case, 
this is done in the next section. 

We conclude this section with mentioning an approach to finding solutions to the bifurcation 
equation (|8.6|) for any n. For a fixed n, we define the first reduced energy E\(ip) := £\(^),a), 
where a = A% + a, with A%(x) := %Jx and a(ip) = {M + \4>\ 2 )- 1 Im(^V^n^) (see © and 
Critical points of this energy solve the equation F(X, ip) = 0. 

Next, we introduce the finite dimensional effective Ginzburg-Landau energy 

e\(v) := E\(v + w(X,v)). 

It is a straightforward to show that 

(i) e\(v) has a critical point vq iff E\(u) has a critical point uq = vq + w(X,vq); 

(ii) Critical points, vo, of e\(v) solve the equation (|8.6|) ; 

(iii) e\(v) is gauge invariant, e\(e lS v) = e\(v). 

One can use e\{v) to investigate solutions of the equation (|8.6p for any n. 
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9 Proof of Theorem 11.11 

In this section we look at the case n = 1, and look for solutions near the trivial solution. For 
convenience we drop the (super)index n = 1 from the notation. Recall that ip is a non-zero 
element in the nullspace of the operator L n — n acting on M* n (r) . Since by Proposition 16.11 
this nullspace is a one-dimensional complex subspace for n = 1, the Abrikosov function, f3(ip ), 
defined in (|5.4p . depends only on r. Therefore we write (3{t) = (3(ipo), so that 

Of S . (I^OI 4 ) /Q -1 \ 

^ = (Wf (9 - 1} 

We begin with the following result which gives the existence and uniqueness of the Abrikosov 
lattices. 

Theorem 9.1. For every t there exist e > and a branch, (A s , ip s , a s ), s € [0, y/e), of nontrivial 
solutions of the reseated Ginzburg- Landau equations (|4.3[) . unique modulo the global gauge sym- 
metry ( apart from the trivial solution (1, 0, Aq) ) in a sufficiently small neighbourhood of (1, 0, Aq) 
m R x J^(t) x J^(t), and such that 

'X s = 1+ 5 a(s 2 ), 

ip a = stPo + sg^s 2 ), (9.2) 
a s = A + g a (s 2 ), 

where (L — l)"0o = 0, gy, is orthogonal to null(L — 1), g\ : [0, e) — > R ; : [0, e) — > ^°(t), and 
g a : [0, e) — ¥ J4?(t) are real-analytic functions such that g\(0) = 0, <?i/>(0) = 0, g a {0) = and 



SX(0) = 



(l^ol 2 )- (9.3) 



Proof. The proof of this theorem is a slight modification of a standard result from bifurcation 
theory. Our goal is to solve the equation (|8.6[) for A. Since the projection P, defined there, is 
rank one and self-adjoint, we have 

P ^ = ITTH (^0,^)^0, with G null^F(Ao,0). (9.4) 

We can therefore view the function 7 in the bifurcation equation (18.6[) as a map 7 : R x C —> C, 
where 

7 (A, s) = (V'o, F(A, s^o + w(X, s^ )). (9.5) 

We now show that 7(A, s) G R. Since the projection Q is self-adjoint, Qw(X, v) — w(X, v), w(X, v) 
solves QF(X, v + w) = and v = sifto, we have 

(w(X, sip ),F(X, sipo + w(X, sipa))) = (iu(A, sip ), QF(X, sip + w(X, sipa))) = 0. 

Therefore, for s ^ 0, 

(ip ,F(\,sip + ®(\,sipa))) = s-^sipo + w(X,sip ),F(X,sipo + w(X,sip ))}, 

and this is real by property (jej) of Proposition 17.41 Thus, since by Lemma 18.11 7(A, s) — 
e jargs 7(A, |s|), it therefore suffices to solve the equation 

7o(A,s)=0 (9.6) 
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for the restriction 70 : R x R -> M of the function 7 to R x R, i.e. for real s. Since by 
w(X, sipo) = 0(s 2 ) and therefore (|9.6[) has the trivial branch of solutions s = for all A. Hence 
we factorize 7o(A, s) as 7o(A, s) = S7i(A, s), i.e. we define the function 

71 (A, s) := s- 1 ^, s) = (^ ,F(X, ip + s~ 1 w(X, sip )) (9.7) 

and solve the equation 71 (A, s) — 0. The definition of the function 71 (A, s) implies that it has the 
following properties: 7i(A, s) is real-analytic, 7i(A, — s) = — 7i(A, s), 71(1,0) = and, by (|7.5p 
and (|8.5p . (9a7i(1,0) = — ||^o|| 2 7^ 0. Hence by a standard application of the Implicit Function 
Theorem, there is e > and a real-analytic function <p\ : (— e, e) — > R such that <p\(0) = Xq and 
7(A, s) = with |s| < e if and only if either s = or A = 4>\(s). 

We also note that by the symmetry, 0a(~ s) — 4>\(\s\) = 4>\(s), so 4>\ is an even real-analytic 
function, and therefore must in fact be a function solely of \s\ 2 . We therefore set (j)\{s) = (p\(\/s), 
and so <p\ is real-analytic. 

We now define <f>^, : (— e, e) — > R to be 

W) = \ n n ( 9 - 8 ) 

10 s = 0, 

0^ is also real-analytic and satisfies s4>^{s 2 ) = w{(f)\{s 2 ), si/jo) for any s £ [0, y/e). 

Now, we know that there is a neighbourhood of (Ao,0) in R x null 9^F(Ao, 0) such that in 
that neighbourhood F (A, ip) = if and only if 7(A, s) = where Py = sipo. By taking a smaller 
neighbourhood if necessary, we have proven that P(A, ip) — in that neighbourhood if and only 
if either s = or A = <j)\(s 2 ). If s = 0, we have ip = sipo + scf>^(s 2 ) = which gives the trivial 
solution. In the other case, ip = sipo + scp^{s 2 ). 

The above gives us a neighbourhood of (1,0) in R x J4?(t) such that the only non-trivial 
solutions of the equation (|7.4|l are given by the first two equations in ()9.2j) . We now define 
the function g a (s) = a(ip s ), where, recall, a(ip) is defined in (|7.3|) . This function is real-analytic 
and satisfies g a (— s) = a(—ip s ) — g a (s), and therefore is really a function of s 2 , g a (s 2 )- Define 
a s = Aq + g a (s 2 ). Then (X a ,ip s ,a a ), s € [0,y/e), solve the rescaled Ginzburg-Landau equations 
(Oil . 

We identify (|9.2p with (|5 . 1 [) of Proposition 15.11 with e = s and n = 1. Then (|5.3I) implies 
dO l. □ 

Note that the definition A = (n = 1), the first equation (|9.2I) and the relation (|9.3p imply 
that for (k 2 — \)P{t) + \ > 0, the bifurcated solution exists for b < n 2 , and for (k 2 — i)j8(r)+ h < 
0, it exists for b > n 2 . Thus Theorem 19. 1[ after rescaling to the original variables, implies (I) - 
(II) of Theorem O 

Recall that b — Since the function A s = 1 +g\(s 2 ) given in Theorem 19. 1[ obeys g\(0) = 
and g' x (0) ^ 0, provided (k 2 - |)j9(r) + |^0, the function b = k 2 Aj 1 can be inverted to obtain 
s = s(b). We can define the family (ip s (b)> a s(b) i ^s(6)) of £ T -periodic solutions of the Ginzburg- 
Landau equations parameterized by average magnetic flux b. Since s(b) is real - analytic in b, 
so are ip s (b)i a s(b); A s (&)- This proves (III) of Theorem ll.il 

Due to (|5.13j) . we can express the bifurcation parameter s 2 in terms of b as 

Furthermore, the equations (|5. 15|) implies that the energy of the state (ip s (b), a-s(b) , A s (6)), 

P 6 (r) :=^A(6)(^a(6) 5 a«(6))) ( 9 - 10 ) 



AbrikosovLattices, December 4, 2011 



15 



where we display the dependence on r, coming through the solution (ip Sl a s ), has the following 
form 

The next result addresses the nature of dependence of solution on r. 
Lemma 9.2. (X Sl ?p s ,a s ) depend smoothly on t± := Rer and r 2 := Imr. 

Proof. By the above we can write f/'s = s^o + w(A s ,s) and a s — Aq + a(w(X s , s), s), where 
a(w, s) is given by 

a(w,s) = (M + |sVo + w] 2 )' 1 Im((s^ + w)V Ao (s^ + w)), (9.12) 

w(X, s) solves the the equation 

(L-X)w = -Q[2ia(w,s)-V Ao (sV^ + w) + |a(w,s)| 2 (# + ^) + K 2 |(sV'o + w)| 2 (s^o + w)]^ (9.13) 

and A s solves the equation j(X,s) — (see (|9.6|) ). Here L is the restriction of L to RanQ. 
From their explicit expressions we see that ipo, P and Q are smooth in r,. Differentiating (19. 13)) 
with respect to n and solving the resulting linear equation for d n w s , it is not hard to convince 
oneself that w s is differentiable in Tj. Repeating this proceedure, one sees that w s is smooth in 
r,-. Therefore ip s and a s are smooth in r,-. For the same reason, 7(A, s) := PF(X, sipo + w(X, s)) 
is smooth in 73 and therefore so is A s . □ 

Corollary 9.3. Eb(r) depends smoothly on t\ := Rer and r 2 '■— Imr. 

Proof. To get rid of dependence of the domain of integration in (|4. 2[) on r, we change the varables 
of integration in (j4.2[) as x — m T y, where m T = (Vlmr) -1 f J Imr 1' wn i cn reduces the 
integral to the one over the unit square cell. Since the rescaled functions, 

{$s(x) = ip s (m T x), ^ 
1 a s (x) = m l T a s {m T x), 

defined on a r-independent square lattice, are still smooth in t\ :— Rer and r 2 := Imr, the 
result follows. □ 

The next result establishes a relation between the minimizers of the energy and Abrikosov 
function. 

Theorem 9.4. In the case K > -^=, the minimizers, Tb, of r i— > J5&(t) are related to the 

minimizer, r*, o/ /3(t), as r^ — r* = 0{p}^ 2 ), where /i := k 2 — 6. In particular, Tb — > r* as 
6^k 2 . 

Proof. To prove the theorem we note that Eb(r) is of the form Eb(r) — eo+ei/i+e 2 (r)/i 2 +0(/i 3 ). 
The first two terms are constant in r, so we consider Eb(r) = e 2 (r) + O(fi). Tb is also the 
minimizer of r H> Eb(r) and r* , of e2 (r ) . We have the expansions -E& (r* ) — Eb(rb) = \E h (r^ ) (r* — 

r 6 ) 2 + 0((r» - r 6 ) 3 ) and E 6 (r,) - £ b (r h ) = - |e' 2 '(T&)(T* - r b ) 2 + 0((r, - r h ) 3 ) + 0(/x), which 
imply the desired result. □ 

The following result was discovered numerically in the physics literature and proven in [2] 
using earlier result of [23] : 
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Theorem 9.5. The function /3(r) has exactly two critical points, r = e ,lr / 3 and r = e l,r / 2 . The 
first is minimum, whereas the second is a maximum. 

Theorems I9.4I and I9.5I imply the remaining statement, (IV), of Theorem ll.il □ 
Remarks. 1) The applied magnetic field is given by ho = ^dbEtir). From (|9.1ip we have 
formally 

= t+Y(l*o| 3 ) + 0((K 2 -i») 2 ). (9.16) 

2) The proof of Theorem 19.11 gives in fact the following abstract result. 

Theorem 9.6. Let X and Y be complex Hilbert spaces, with X a dense subset ofY , and consider 
a map F : M x X — > Y that is analytic as a map between real Banach spaces. Suppose that for 
some Ao € K, the following conditions are satisfied: 

(1 ) F(X, 0) = for all A € M, 

(2) d^F(\o,0) is self-adjoint and has an isolated eigenvalue at of (geometric) multiplicity 1, 

(3) For non-zero tp € null^F(A , 0), (ifo, d x ^F(X , 0)Vo) ^ 0, 

(4) For all ael, F(X, e ia ip) = e ia F(X,tp). 

(5) For all tp G X, (tp, F(\, tp)) G K. 

Then (Ao, 0) is a bifurcation point of the equation F(X, ifi) — 0, in the sense that there is a family 
of non-trivial solutions, (A s , ip s ), for s € [0, y/e), unique modulo the global gauge symmetry (apart 
from the trivial solution (1,0)) in a neighbourhood of (Ao,0) in K x X. Moreover, this family 
has the form 

|V = s ^o + s<j)^(s 2 ). 

Here ipQ € null 9^F(Ao, 0), and 4>\ : [0, e) — > M and : [0,e) — > X are unique real-analytic 
functions, such that 4>\(0) = Ao, (0) = 0. 

A Fixing the Gauge 

We provide here an alternate proof of Proposition 14.11 largely based on ideas in [T5] . We begin 
by defining the function B : M — >• R to be 



B(C) 



- f coil A(Z,Qd£. 
r Jo 



It is clear that b — J Q TT2 B(Q d(. A calculation shows that £?(£ + rr 2 ) = B(Q. 
We now define P = (Pi, P 2 ) : K 2 -> K 2 to be 



Pi(a:) = 6x 2 - / * B(()dC, 



o 



ft (a) 



/ curl A(£, x 2 ) d£H B(x 2 ). 
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A calculation shows that P is doubly-periodic with respect to C. 
We now define 77' : K 2 -> M to be 

t/(x) = Jxix 2 - / * Ai(Z,0)dt- [ * A 2 (xi,()- P 2 (x 1 ,()d(. 
1 Jo Jo 

7/ satisfies Vrf = - A + A Q + P. 

Now let 77" be a doubly-periodic solution of the equation At/' = — div P. Also let C = 
(Ci,C 2 ) be given by 

C =-H/n (P + V '" )<te ' 

where f2 is any fundamental cell, and set rj'" — C\X\ + C 2 x 2 - 

We claim that 77 = 77' + 77" + 77"' is such that A + V77 satisfies (i) - (iii) of the proposition. 
We first note that A + Vr] = A- A + A Q + P + V77" + C. By the above, A' = P + V77" + C is 
periodic. We also calculate that div A' = divP + A77" = 0. Finally J Q A' = f n (P + Vr)-C) = 0. 

All that remains is to prove (iv). This will follow from a gauge transformation and translation 
of the state. We note that 

Ao(x + t) + A'(x + t) = Ao(x) + A'(x) + ^ ( ^ 

This means that A (x + t)+ A'(x + 1) = A (x) + A'(x) + Vg t {x), where g t (x) = \t A x + C t for 
some constant C t - To establish (iv), we need to have it so that C t = for t = r, rr. First let I 
be such that r A I — — ^ and rr A I — — . This I exists as it is the solution to the matrix 
equation 



r \ ( h 



-rr-, rr, : \ I ■> ! I 



and the determinant of the matrix is just r 2 r 2 , which is non-zero because (r, 0) and rr form 
a basis of the lattice. Let ((x) = |Z A x. A straight forward calculation then shows that 
e l C( x )i[)(x + I) satisfies (iv) and that A(x + I) + V((x) still satisfies (i) through (iii). This proves 
the proposition. 
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